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ABSTRACT 

Two dimensional calculations of the evolution of remnants generated by the 
strong mechanical energy deposited by stellar clusters in dwarf galaxies (M ~ 
10^ — 10^° Mq) are presented. The evolution is followed for times longer than 
both the blowout time and the presumed span of energy injection generated by 
a coeval massive stellar cluster. The remnants are shown to end up wrapping 
around the central region of the host galaxy, while growing to kpc-scale 
dimensions. Properties of the remnants such as luminosity, size, swept up mass, 
and expansion speed are given as a function of time for all calculated cases. 

The final fate of the swept-up galactic gas and of the matter processed 
by the central starburst is shown to be highly-dependent on the properties of 
the low density galactic halo. Superbubbles powered by star clusters, with 
properties similar to those inferred from the observations, slow down in the 
presence of an extended halo to expansion speeds smaller than the host galaxy 
escape velocity. Values of the critical luminosity required for the superbubbles 
to reach the edge of the galaxies with a speed comparable to the escape speed 
are derived analytically and numerically. The critical luminosities are larger 
than those in the detected sources and thus, the superbubbles in amorphous 
dwarf galaxies must have already undergone blowout and are presently evolving 
into an extended low density halo. This will inhibit the loss of the swept-up and 
processed matter from the galaxy. 



-3- 
1. Introduction 

Massive stellar clusters violently eject their processed matter through supernova 
(SN) explosions and strong stellar winds, leading to the formation of large-scale remnants 
commonly known as superbubbles. These, when generated in dwarf galaxies, are thought 
to acquire a dimension that rapidly exceeds, at least along the galaxy minor axis, the 
scale height of the mass distribution. When this occurs, the remnants are believed to blow 
out of their host galaxy i.e. as the shell of swept up interstellar matter (ISM) accelerates 
outwards and fragments via Rayleigh-Taylor (R-T) instabilities, it allows the processed 
matter to freely stream as a wind into the intergalactic medium (IGM). Furthermore, the 
most energetic events are thought to cause not only the loss of their newly processed metals 
but to be able to strip dwarf galaxies of their entire ISM (see De Young & Heckman 1994 
and references therein). 

Here however, we show that despite the fact that remnants easely manage to blowout 
or burst out of the disk matter distribution, under a wide variety of circumstances the loss 
of matter could be inhibited by the hydrodynamic response of the rapidly streaming wind 
gas to the conditions found in the galaxy halo. 

Analytical and numerical hydrodynamic studies of the remnants caused by a 
strong stellar energy input in a plane-stratified atmosphere have analyzed the blowout 
phenomenon (Kompaneets 1960; Tenorio-Tagle & Bodenheimer 1988; Mac Low et al. 
1989; Igumentshehev, Tutukov & Shustov 1990; Tenorio-Tagle, Rozyczka & Bodenheimer 
1990). For a recent account see Koo & McKee (1992) who derived the characteristic wind 
luminosity (L^) that determines whether a bubble remnant is able to blowout of the disk 
configuration. Lh is the luminosity required for the expanding shell of swept up matter to 
have a sonic velocity when it reaches a disk scale height. Lb — 2.4 x 10^^P4iffoo^s,io ©rg s"^, 
where P4 = Po/(10'^/c) cm~'^K, i^ioo is the disk scale height in units of 100 pc and ag^io is 
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the disk sound speed in units of 10 kms~^. Koo & McKee (1992) estimated also that if the 
cluster luminosity {L^) exceeds SL^ the bubble would expand supersonically when reaching a 
scale height and it will begin to accelerate to blowout of the disk into the halo of the galaxy. 
Once the bubble blows out, there are two possible end solutions depending on the wind 
luminosity. For » blowout would cause the venting of the hot bubble interior out of 
the planar atmosphere. Lower wind luminosities (3 < L^,/Lft < 10) result in the formation 
of stable hydrodynamical jets as the processed matter is re-shocked and accelerated to 
supersonic speeds through a nozzle generated by the shocked disk matter. In this particular 
case there is no doubt the processed matter is due to escape the galactic system. The more 
energetic events ( >> ) destroy upon expansion the possibility of a nozzle formation 
and the hot gas in such cases acquires a maximum speed similar to its own velocity of 
sound, as it streams through the fragmented shell into the halo. In dwarf galaxies, these 
remnants have recently been investigated by several groups. In particular, amorphous 
dwarf galaxies have been selected as ideal laboratories to study the impact of starbursts 
on the ISM, and to investigate the possible ejection and loss of their ISM and processed 
matter into the IGM. These galaxies are ideal laboratories because they have an ISM mass 
in the range 10^ - 10^ M© and a concentrated central starburst, as opposed to Magellanic 
Irregulars where the star formation seems to extend over the whole body of the galaxies 
(Gallagher & Hunter 1987; Marlowe, Heckman, Wise & Schommer 1995 and references 
therein). Theoretically these issues have been address by De Young & Gallagher (1990) 
who estimated that about 2/3 of the processed matter ejected by supernova will escape the 
galaxy. This conclusion was reached from an analysis of 2-D numerical calculations in a 1 
kpc^ grid, after comparing the processed matter expansion speed with the escape velocity 
of their selected galaxy. In a further study, De Young & Heckman (1994), who considered 
dwarf galaxies with masses in the range lO'^ — 10^ M© and physical sizes of 1 - 3 kpc, found 
that the degree of flatness of the assumed density distribution is an important parameter. 
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This is because after shock acceleration and through the blowout phenomena, the energy 
of the starburst tends to escape along the galaxy minor axis imparting little momentum to 
the gas in the outer portions of the disk matter distributions. Thick disks however, were 
found to be completely disrupted and accelerated to speeds larger than the escape velocity 
of galaxies. Note however, that the approach of De Young & Heckman does not include 
the contribution of dark matter to the gravitational potential of the dwarf galaxy, nor the 
existense of an extended halo component. Note also that the range of energies used spans 
to up to 10^^ erg s~^, although the maximum assumed energy was never allowed to surpass 
the number of SN expected if the whole galaxy ISM (Mjsm) was turned into stars following 
a Salpeter IMF i.e. the maximum used number of supernovae N^jv = M/^m/IOO. It is 
important to note that the best studied dwarf galaxies have an estimated mechanical energy 
deposition in the range of 10^° - 10^^ erg s~^ (Marlowe et al. 1995). This imphes a SN rate 
of 3 X (10~*^ — 10"^) per year and a total expected energy injection of 1.3 x (10^^ — 10^^) 
erg if the energy input rate is assumed to last for 4 x 10^ yr, time after which all massive 
stars from a coeval starburst would have exploded as SN. 

Unfortunately, there is little information about the gaseous structure of dwarf galaxies. 
However, in all well studied cases the size of the region photoionized by the cluster of 
massive stars, exceeds the size of the superbubbles present in these galaxies (see Marlowe 
et al. 1995; their tables 4 and 5). Furthermore, the few that have been looked at in HI 
(see for example Hoffman et al. 1993; Hunter, Van Worden and Gallagher, 1994; Meurer 
1994) display large envelopes, several times larger than the optical size of the galaxies. 
One can obviously argue that it may just be a projection effect, however, the galaxies 
under consideration are of the " amorphous" type and thus it is rather unlikely that all 
observed matter is confined to a narrow disk. It thus seems perhaps more likely that the 
superbubbles observed in these galaxies (Marlowe et al. , 1995; Martin 1996) which present 
a continuous horseshoe shape attached to the central star forming region and which have 
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dimensions ~ 1 kpc and expansion velocities ~ 50 kms ^are presently evolving into the 
extended halos of their host galaxies. 

Here we present two dimensional calculations of multi-supernova remnants evolving in 
dwarf galaxies (see sections 2 and 3) . We find a good agreement with the above mentioned 
studies although our conclusion (see sections 4 and 5) is that the final fate of the swept-up 
and processed matter is highly-dependent on the properties of the previously ignored 
galactic gaseous halos. 

2. Dwarf GalzLxies 

2.1. The ambient gas density distribution 

Massive galaxies (Moai > 10''^° M© ), such as the Milky Way, are capable of retaining a 
hot {Tcai ~ a few lO^K) and extended massive halo completely bound to their gravitational 
potential. This is frequently re-heated by a collection of superbubbles, caused by aging OB 
associations, that blowout of the plane HI disk density stratification. Less massive galaxies, 
on the other hand, are unable to keep such massive and hot halos. Instead, the smaller the 
galaxy total mass, the cooler and less extended their possible resultant halos become. Thus 
these are usually idealized as a low density turbulent medium (or clouds) that extends up 
to a maximum radius Rq at which its characteristic velocity dispersion exceeds the galaxy 
escape speed. As stated by Suchkov et al. (1994, hereafter SBHL), "in a real galaxy, much 
of the hydrostatic support of the interstellar matter is provided by the random motions of 
interstellar clouds" . 

Following Li & Ikeuchi (1992), Tomisaka & Ikeuchi (1988, hereafter TI), Tomisaka & 
Bregman (1993, hereafter TB), and SBHL, we model our dwarf galaxies with three different 



- 7- 



isothermal components. 

Pg = PNM + PlM+ < Phalo > (1) 

where pnm, Pim and < phaio > are the gas densities of the neutral and ionized components 
and the mean density of the extended turbulent halo material, respectively. The model 
corresponds to the bound-cool type halo in the notation of Li & Ikeuchi (1992). The total 
gas pressure is then given by 
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where Cnm, Cjm and Chaio are the sound velocities of the neutral and ionized gas and the 
velocity dispersion of the extended halo component. 7 is the ratio of specific heats. 

Let Rb he the maximum extent of the stars and dark matter, and Vb, the escape 
velocity at this surface. Following TI, and TB we assume that a fraction (1 — e^) of the 
radial component of gravity is balanced by pressure gradients. Then the initial gas density 
distribution follows from the relation: 
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where po is the total gas density at the galactic center, = Pnm,o/po, c^im = Pim,o/po 
and afiaio =< Phaiofi > / Po are the density ratios of the various components to the total 
density at the galactic center. The function x equals 



X = FH - e^F{r) - (1 - e')F{0), 



(4) 



where r = y/x'^ + y'^ and 00 = \/ x'^ + y'^ + z'^ are the cylindrical and spherical radii, and the 
function F is defined by equation (|10D- The rotation velocity of the ISM component comes 
from the relation: 

^ = e'{9l + 9ir/\ (5) 
where Qx and Qy are the x and y components of the gravity field. 
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We include a transition from a rotating disk component to a non-rotating spherical 
halo and radial gradients in the disk density distribution by assuming 



where Hr and Hz are the characteristic scale heights along and perpendicular to the plane 
of the galaxy. Then @ is not self-consistent if a hydrostatic gaseous disk is assumed, 
but it is asymptotically correct for a hydrostatic gaseous halo (TI, SBHL). Typical dwarf 
galaxy density distributions assuming a total mass Mb = 10^° Mq (models A and B) and 
10^ Mq (Models C) and a 10% gaseous component are shown in the Figure 1 (see also 
Table 1). Note that the last contour drawed represents, in each case, the gaseous edge of 
the galaxy {Rq)- i-e. material further out with a velocity dispersion Chaio would exceed 
the escape velocity and thus it would not be bound to the system. Figure 2 gives a global 
impression of how the galactic ISM is distributed for different values of the assumed velocity 
dispersion of the extended component (also used to identify the models; see Table 1) and 
thus it shows how for larger values of Chaio, the smaller, or more compact that the bound 
gaseous component results. The value of PextiRc)/^ varies in our models within a range 
0.8 cm"3 K (Model B60) - 5.4 x 10"^ cm-^ K (Model A120). This pressure should balance 
the extragalactic gas pressure and in fact covers the range often used in the literature (eg. 
Babul & Rees 1992). 

Several calculations were also made for different ratios of Mjsm/Mb, as well as for 
different dimensions of the dark matter distribution (see section 5). 




(6) 



2.2. The gravitational field 



We have assumed that the main contribution to the gravitational field results from 
spherically symmetric stellar and dark matter components, disregarding the contribution of 
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a stellar disk component, and the self-gravity of the gaseous matter. A King model was 
assumed for the background stars and a distribution inversely proportional to the square of 
the radius has been adopted for a dark matter halo: 

Pi , P2 



2 ' 



(7) 



where Ri and i?2, Pi and p2 are characteristic scales and densities of the star and dark 
matter components, respectively. 

The total mass of these components may be expressed as 
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where xb — Rb/Ri, Mi — AnpiRf, ub — Rb/R2, M2 — A'Kp2R\. We fix a constant at 
the gravity potential by assuming that it becomes Newtonian at the galactic boundary 
Rb- Then the gravity potential corresponding to the adopted total density distribution is 
expressed as follows: 

$(r,^) = -^FH, (9) 
where function F{uj), is defined by the expression (see also SBHL) 
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F{u)^-^, u>Rb, 
where — uj/Ri, — uj/ R2. The value of Vb is given by 
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whereas the escape velocity at the current position (r, z) is 



Vesc = VBjF{iU 



(12) 



The X, y, z components of gravity for u < Rb derived from (H) and (0) are 
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For larger distances uj > Rb they are simply 
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(13) 
(14) 
(15) 
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The Calculations 



Here we present two-dimensional calculations of the evolution of the remnants 
generated by multi-supernova explosions from an aging massive stellar cluster evolving in 
a low-mass galaxy {Mb < 10^'^Mq ). The remnant evolution is calculated by means of 
the thin layer approximation (Kompaneets 1960) a method developed in 2D and 3D by 
Bisnovatyi-Kogan & Blinnikov (1982); Mac Low & McCray (1988); Bisnovatyi-Kogan & 
Silich (1991); Palous (1992); Silich (1992) and Silich et al. (1996), and now applied to a 
variety of cases (Bisnovatyi-Kogan & Silich 1995). The evolution is followed up to 60 - 
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170 Myr, long past the time at which the SN activity from the evolved - assumed coeval - 
massive cluster has come to an end. 



3.1. The energy input 

Based on the observations of dwarf galaxies (see Marlowe et al. 1995 and references 
therein) and the synthetic properties of starburst galaxies derived by Leitherer & Heckman 
(1995), the energy input rates for the calculations were set within the range = 10^° - 10^^ 
erg s~^, although in section 5 much more energetic clusters are considered. The assumed 
constant energy input rate depends on the total energy released by massive stars {Eburst) 
and the duration of the supernova phase (hurst) and can be written as 

Lss = 3.17 X ^Q^ Eburst/lO^'ergs ^ 

Here the energy supply rate is L38 — / 10^* erg s~^. We have also assumed that massive 
stars release an average msN — lOM© during their explosion, and that this is dispersed into 
the hot superbubble interior. Thus the rate of mass ejecta due to SN explosions, Msn-, is 
assumed to be 

■ NsNmsN L^msN .r,n\ 
Msn = — = — ^; i^Uj 

tburst J^SN 

where Esn is the energy released by each SN explosion (typically Esn = lO^^erg). 



3.2. Initial conditions 

All the calculations start with a small adiabatic spherical bubble of initial radius Re at 
the time 
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75(7 - 1) 



(21) 
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The initial expansion velocity, thermal energy, and mass within the bubble are 

uo = 0.6^, (22) 
14 

Eth(to) = - — —LJo, (23) 
97 + 5 

Minito) = MsNto- (24) 

Table 1 provides the main parameters of the simulated galaxies, as well as the energetics 
of their assumed bursts of star formation. Several parameters were unchanged in the 
calculations e.g. the metallicity of the galactic gas has been assumed to be ^ = Z/Zq = 0.3; 
the duration of the supernova phase tb^rst = 4:0 Myr; the constant Cq = 0.9 in (eq H). The 
parameter n^Mfi = 20 cm~^, (see (^)) the characteristic scale heights for the gas density 
distribution have been taken as = Ikpc, = 0.5Rb, whereas the characteristic scale 
height for the King's mass distribution was taken to be Ri = 0.75 kpc, and the parameter 
i?2 in equation (^ has been kept as half the truncation radius Rb- The temperature of the 
ionized component was taken as 2 x 10^ throughout the calculations. To satisfy the high 
pressure condition at the starburst region we also set the neutral gas temperature equals to 
6 X 10^ K. We relaxed the high pressure constraint for the low mass galaxies (models C), 
for which a temperature for the neutral gas of 100 K was assumed. 

The method of solution and the procedure to calculate various remnant properties such 
as the bolometric and X-ray luminosities are given in Appendices A and B. 



4. A typical evolutionary sequence 

Figure 3 shows the evolutionary sequence of models AlOO (see Table 1) for which 
the initial ISM density follows the distribution shown in Figure lb with a central density 
value of 20.2 cm~^and a halo velocity dispersion of 100 kms~^. The central starburst was 
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assumed to deposit 7.9 x lO^'^ erg during the first 40 Myr of evolution. As expected, the 
remnant grows at first elongated in the direction perpendicular to the disk of the galaxy 
and blows out into the halo. Soon however, as a result of blowout a secondary shell of 
newly swept up halo gas forms. This results from the stream of the high pressure processed 
matter that bursts out of the superbubble interior while causing a symmetrical and larger 
twin superbubble(see Figure 3a). This immediately begins to wrap around and ends up 
surrounding the inner densest part of the disk as it grows to large dimensions into the halo. 
Note that more elongated structures result if steeper density distributions are assumed (see 
Figure 3b). At later times, the leading shock manages to plough through the outer less 
dense sectors of the central disk and merges with its symmetrical counterpart to form a 
continuous structure and a more spherical single giant superbubble. During the process of 
shock merging, a considerable amount of the shell mass is left behind strangling a dense 
toroid carved in the initial density distribution. All of this matter stops then participating 
in the general outward motion and is instead strongly compressed towards the galaxy plane 
to remain completely engulfed by the almost spherical giant superbubble. 

Figure 4 shows the total amount of matter swept by the remnant as a function of time. 
This steadily grows to values of several times ~ lO^Mo with most of the material being 
collected as the remnant evolves into the halo. The sharp decrease shown in Figure 4a at 
about 25 xlO^ yr is due to the process of outer shock merging above described. Figure 4 
also shows the total amount of matter evaporated from the outer shell into the remnant 
interior, and the amount of material injected by the collection of SN explosions. In all 
cases, the evaporated mass amounts to only a few percent of the swept up mass but it 
constitutes the most significant mass input into the superbubble interior; at least an order 
of magnitude larger than that provided by the sequential SNe. 

Meanwhile, the hot bubble interior, as well as the shell of swept up matter, cause the 
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X-ray emission from the nebulae. Figure 4e displays the total X-ray emission (see Appendix 
B), and the shell contributions to this luminosity. The bolometric bubble luminosity is 
shown in Figure 4d. In the conditions of the A 100 model the remnant evolves mainly in the 
radiative mode and thus the shell contribution to the X-ray emission is negligible compared 
to the X-ray emission from the hot superbubble interior. 

5. The fate of the processed matter 

The central issue however, is whether or not one can estimate without uncertainties the 
possibility of mass ejection from the galaxy. For this, the obvious and clear cut procedure is 
to compare the fastest expansion speed of the remnant (say along the symmetry axis) with 
the galaxy escape velocity at different radii, or evolutionary times. The maximum speed 
(Ug-rp) of the sector of the remnant evolving into the steepest density gradient shows at first 
a continuous deceleration (see Figure 5), for i?'s less than 1 kpc. Then during blowout 
into the galaxy halo, it suffers a sudden increase towards large speeds surpasing even the 
galaxy escape velocity values (dashed lines in Figure 5). The dynamic halo however, makes 
eventually an impact on the evolving superbubble forcing it to slow down again, once a 
significant amount of halo material has been incorporated into the remnant shell. These 
trends are shown in Figure 5 where the largest expansion remnant speeds (solid lines) for 
all A and B cases (see Table 1, and Figures 1 and 2) are compared with the galaxy escape 
velocity. Note that the maximum expansion speed after blowout into the halo, as well as the 
slowest deceleration occurs for the more spread out distributions of matter or less compact 
galaxies considered; which are those with the smallest Chaio values. All superbubbles attain 
large expansion speeds, even larger than Vesc, however, they all end up with a velocity well 
below the galaxy escape speed. 

An important point to notice here is that just as old supernova remnants are 
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disrupted once their expansion speed drops below the ISM random speed of motions, 
once a superbubble expansion becomes slower than the assumed velocity dispersion of 
the host galaxy halo, it should begin to be disrupted and eventually loose its identity. 
Giant, kpc-scale, filaments associated with the host galaxy, may be the only long lasting 
trace of the energetics from old stellar generations and their disrupted giant remnants. 
Thus important during and after disruption is to decide how much further can the hot 
superbubble interior expand. Note that its expansion would lead to further cooling and this 
to a further drop in pressure. The final fate, i. e. whether it will be lost or trapped by the 
galaxy, can be addressed by comparing its pressure {Pm) to the host galaxy outer boundary 
pressure Pq = Pexti^a)- With this aim in mind, the calculations were continued for times 
that largely exceeded the onset of shell disruption (ie. the time when the expansion speed 
becomes smaller than the random speed of motions in the halo), and which provided us 
with a handling on the run of Pj„. This, as shown in Figure 6, in most cases becomes 
smaller than Pq with the implication that the processed matter will remain bound to the 
system. For the case A60, although Pj„ exceeds Pq, it does not exceed the local pressure 
value at the shell stall radius (dotted line in Figure 6) , and this leads to the contraction of 
the superbubble and the retention of the processed material. 

Table 1 indicates several other calculated cases for which we assumed for example 
a different proportion of dark matter with respect to Mjsm (D and E cases), or runs 
without thermal evaporation (A* 100), or with a lower total mass (C cases). The results are 
nevertheless very similar to the above described cases. 

Clearly, the expulsion of matter from a dwarf galaxy will only be possible if one could 
freely increase the energy input rate (L). The problem however, is to know which is the 
correct or sufficient amount for each of the galaxies. A good formalism could be device if 
one notices from Figures 1 and 2 that the extended low density halos are in fact rather 
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uniformly spread out and that they constitute a large fraction of the galaxies ISM. Thus, 
given the continuous energy input rate (at least for the first 40 Myr) one would expect that 
the relationships developed for the evolution of stellar wind bubbles in a uniform density 
medium (see Weaver et al. 1977) will give reasonable estimates. The functions R(t) and 
dR/dt lead to the time independent expression 

R = i,.e(hif("y\p,. (25) 



no J \kms~ 

iFiom this we would like to know the amount of energy required for a superbubble to 
reach the outer radius of a galaxy with a speed comparable to the escape velocity of the 
system. We shall call this the critical luminosity {Lent)- Thus, if one sets R = Rg, u = Y^sc 
and the mean number density no ~ < rihaio > = IMism/^R^I^h] where f is the fraction 
of the galaxy ISM swept during the evolution of the superbubble and [Ih is the mean mass 
per particle, it leads to: 



,_.,T5xlO-(/^)(^)7|i)" (26) 



where Lcwt.ss = LcHt/l^ erg. 



A similar expression can also be found for late evolutionary times (t > iburst) for which 
an instantaneous energy release may seem a better approximation. Thus, from Chevalier 
(1974) and Blinnikov et al. (1982) one obtains: 

R = 22('-^)"'\r PC, (27) 

where E^q is the total amount of energy deposited by the star cluster in units of 10^° erg, 
and ts is the evolutionary time in 10^ years. For the fraction of the total input energy which 
is transformed into the remnant thermal energy e = 5/11 (Weaver et al. 1977). Relation 
(p7|) leads to the critical mechanical luminosity Lcrit = Ehurst/tburst'- 
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Both relationships when the fraction of the swept up interstellar mass f was set to 0.7 are 
displayed in Figure 7a,b for galaxies with an Mjsm = 10^ Mq and 10*^ Mq , respectively. 

Figure 7a (solid line) shows the Lcrit trend for = 5 kpc and the dashed line for 
= 10 kpc (cases A and B). The critical luminosity, although by definition it is sufficient to 
cause the ejection of processed matter out of the galaxy, it may only do so, depending on 
the host galaxy mass distribution, after a long evolutionary time. In some cases this can be 
larger than thirst and in fact much larger than the expected HIT region lifetime (~ 10^ yr). 
Several of our numerical calculations in search of Lcrit are also shown in Figure 7. Diamonds 
in Figure 7a represent the numerical results for galaxies with R^ = 5 kpc, and triangles for 
Rb = 10 kpc models. In all cases the time required for the superbubble to reach the edge 
of the galaxy are indicated. 

There is a better than 30% agreement between the derived formulae and the numerical 
results. For superbubbles that require a time smaller or comparable to the starburst time 
hurst to reach the galaxy edge, our numerical estimates agree with relation (^), and for 
those that take a time longer than Uurst equation (p8|) leads to an asymptotic value of the 
critical energy input rate Lcrit- 

An extended halo can therefore inhibit the loss of the swept-up and processed matter 
from a galaxy even if the remnants undergo blowout from the dense central gaseous 
component into the halo. Note also that the more compact ISM distributions demand more 
powerful stellar clusters (a larger Lcrit) to allow their superbubbles to reach the outskirts of 
galaxies and eject their metals into the IGM. One can therefore conclude that the presence 
of an extended halo may dominate the evolution of the remnants, and that the final fate of 
the swept-up galactic gas and of the matter processed by the starburst, basically depend on 
the properties of the low density galactic halo rather than on the parameters of the gaseous 
central disk. Important is also the extent of the dark matter component: the more compact 
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this is, the larger the energy required to reach the critical luminosity. 

Note also that all calculations with an energy input rate < 10^^ erg s~^ and a halo 
velocity dispersion Chaio > 60 kms~^, as inferred from the observations of amorphous 
dwarf galaxies (Marlowe etal. 1995, Martin 1996, etc.), all lie below the curve of critical 
luminosity, and thus have led to superbubbles expanding with speeds smaller than the 
random speeds of motions of their host galaxy halos, implying their disruption long before 
they reach the edges of their galaxies, and thus have led to no loss of their processed 
matter. Such galaxies in order to loose their metals into the intergalactic medium, require 
of energetic star burst, more energetic than what they presently have. One could in this 
respect be tempted to think of a larger energy, larger than observed, being deposited in an 
extended dwarf galaxy and thus more easely surpase the critical luminosity value. However, 
the time required for such superbubbles to reach the galaxy edge and eject their metals into 
the IGM would exceed by large factors the HIT region hfetime and thus most probably be 
undetectable events. 



6. Conclusions 

Based on the scarce but definite observations of amorphous dwarf galaxies, which show 
that the HI distribution outgrows by large factors the optical size of these galaxies, we have 
constructed realistic models which account for an extended cool halo component bound to 
the system. Models of amorphous galaxies with an Mjsm ~ 10^ Mq and an energy input 
up to 10"^^ erg s~^, as inferred from the observations, lead to blowout from the disk matter 
stratification, allowing the superbubbles to vent their hot matter into the galaxy halo. The 
outfiowing gas immediately searches and follows the steepest density gradients affecting 
a large almost spherical volume. We have shown here that in the case of small systems, 
superbubbles wrap around and envelope the densest central region of the host galaxy. Such 
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remnants are strongly decelerated as they evolve into the halo, sufficiently to remain bound 
to their host galaxy. 

Our calculations confirm the results of De Young and Gallagher (1990), and De Young 
and Heckman (1994) during the early superbubble evolution. However, our ultimate 
conclusions are different as the final fate of the remnants basically depends on the properties 
of the galactic halo and on the total energy of the star burst. The presence of the dark matter 
component also enhances the ultimate retention of the processed material. Under these 
circumstances total retention (i.e. expansion speeds < escape velocity) of the remnants and 
thus of the matter processed by the exploding massive stars is unavoidable for galaxies with 
an ISM mass > than a few 10^ and an energy input rate of up to 10^^ erg s~^. 

We have also derived analytical expressions for the critical luminosity required for a 
superbubble to be able to reach the outskirts of a galaxy and expulse the processed matter 
into IGM. The critical limit has been cross-checked with several numerical calculations for 
a large range of galaxy masses, dimensions and thus surface escape velocities. 

^Prom the calculations, we then conclude that the thoroughly studied dwarf galaxies 
with a gaseous mass ~ 10^ Mq and presenting giant remnants produced by the mechanical 
energy (L < 10^^ erg s~^) from massive stars (see Marlowe et al. , 1995; Martin 1996; and 
Meurer 1994) are most hkely to retain their ISM and their processed matter. The observed 
remnants present dimensions smaller than a few kpc and expansion speeds ~ 50 kms~^and 
thus based on our calculations, we predict that they have already undergone blowout 
from their disk matter distribution and that presently are evolving into the extended low 
density halo of their parent galaxy. Our results also predict that these superbubbles will be 
ultimately dispersed by the turbulent motions in the extended galactic halo. 

The consequences of these conclusions are important in a variety of fields and we shall 
analyse these in a forthcomming communication. Clearly, the kind of dwarf galaxies that 



-20- 



we have analysed, are not responsible for the contamination of the IGM. Dwarf galaxies 
with an ISM mass of the order of 10^ Mq will retain their metals unless they undergo an 
overwhelming burst of star formation much larger than those presently observed. We thus 
have to understand why do they present low metalicities. The secret must probably lies on 
the physics of mixing (see Tenorio-Tagle 1996). 

We thank our anonymous referee for many suggestions that allow us to revise and 
improve our knowledge of dwarf galaxies. SAS acknowledges support from the Royal 
Society grant for joint projects with the former Soviet Union, and the RGO and the loA in 
Cambridge for partial financial support. GTT gratefully acknowledges the BBV foundation 
for a visiting professorship at the University of Cambridge. We also thank the hospitahty 
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Appendix A: Main equations 

Following TI, Li & Ikeuclii (1992), TB, SBHL our approach includes three co-existing 
ISM phases: a neutral, an ionized component, and a low density extended halo characterized 
by a turbulent velocity; also thought by other authors as a collection of low densities clouds 
with a velocity dispersion Chaio- 

Mass conservation. The expanding shell caused by the star cluster energy deposition 
accumulates interstellar mass according to 



dt ~ dt dt dt ' 



(Al) 



where the first term represents the mass growth due to swept-up interstellar gas, the second 
one the loss of mass resultant from the development of Ray leigh- Taylor (RT) instabilities 
after the shell blows out of the galaxy disk, and the third one describes possible shell 
self intersections following cusp formation (see below). We have assumed that the swept 
turbulent halo sticks to the shell. Then for any Lagrangian element with a surface area dS 
and a unit vector normal to the shell nj, the terms in equation ([AID t^^e the form: 



d/Xg _ d/iic , d < l^halo > 



dt dt dt 
d/iic 



+ (A2) 



dt 

d < fihalo > 



(PNM + P/M)(Ui - v) • n^dS, Un > 

(A3) 

0, Un < a.. 



dt 



-< Phalo > (U,; - \halo) njdS, (A4) 
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RT 



0, ttsh < 

IJ-L{t=tRT) (A5) 
0, UrT > IJiL{t = tRT) 

— 5{t — t cusp) fJ'L, rem i-^^) 



dt 

An effective sound veiocity = (^nC^M + AC'Im)^^^, wliere An = Pnm/{pnm + Pim), 
A — PimUpnm + Pim)-i is tlie normal component of the shell acceleration, tcusp is the 
time of cusp formation, tar marks the time for the beginning of the acceleration phase, 
Piit = tur) is the mass of the Lagrangian mesh at the beginning of blow-out, and Trt is 
the characteristic growth time for RT instabilities 

where A is the same order of magnitude as the characteristic scale for the ambient gas 
density distribution, A ~ 2Hz (Koo & McKee, 1990). For simplicity we set here A equal to 
the characteristic scale for the gas density distribution at the galactic center A — H^q. 

Momentum conservation. Wc consider four effects leading to changes in the shell 
momentum: the pressure difference between the inside gas and the surrounding medium 
AP, the galactic gravitational force g, the accumulation of the ambient gas momentum, 
and the momentum loss during the RT break-out. We assume that RT clumps separate 
from the shell with the same speed as the current shell expansion velocities. Also that 
intercloud mass accumulation stops wherever the expansion velocity decreases below the 
sound speed a^, and thus the interaction with the ambient intercloud gas modifies only the 
normal component of the momentum but the shell does not accumulate mass any further. 
With these simplifications, the equation of motion for any Lagrangian element with mass 
Hl is (see also Ostriker & McKee 1988, Silich et al. 1996): 

= APndE + + -naio^ - + 
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+I^S - Pic{x, y, z) [(u - v) • iijl^ndE (A8) 
Then the equation for the shell expansion velocity can be written as 

du -lAD , / ^ -l^f^ic , , X -\d\Xhalo , 

= I^L APndS + (v - u)iil + (vcM - nj/i^ + 
+g - pic{x, y, z)ii2^[{u - v) • nJ^ndE, (A9) 

In the current calculations we have assumed that the intercloud gas and halo clouds move 
with the same regular velocities (v = yhaio)- 

The energy balance. The total energy balance for a remnant caused by the strong 
energy input from the central star cluster reads as: 

EsN + Ek,g + Eth,g — Eth + Ek^sh + Eth,sh + Eg + Er^in + Er,sh + Ek,RT: (^H) 

where the left-hand terms are the energy deposited by massive stars (Esn), and the initial 
kinetic and thermal energies of the accumulated ISM {E^^g and Eth,g)- The right-hand 
terms are the bubble interior thermal energy {Efh), the kinetic and the thermal energies of 
the shell {Ek^gh and Eth,sh): the shell gravitational energy (Eg), the energy losses from the 
hot bubble interior (-Br,m) and from the shell {Er,sh), and the kinetic energy of RT clumps 
{Ek,RT)- The derivative of the bubble thermal energy could be expressed as follows 

^ - T -T - p. ^ rA19l 

— ^burst ^in (J/; ' \ ) 

where 

dQ ^ 

^ = E«ndS, (A13) 
i=i 

Lin is the energy loss from the hot bubble interior, and Un is the normal component of the 
expansion velocity. As the main contribution for the interior bubble cooling comes from the 
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outer dense layers, we approximate Lj„ as 

Lin = Lr— = e2LR, (A14) 

where E is the total surface of the remnant, and the index R denotes the radiative part of 
the shell, 62 = S-^/S. Lr is the energy loss from the bubble interior if all the remnant 
becomes radiative (for the calculations of Lr see Silich et al. 1996). 

One can then calculate the bolometric shell luminosity from the total energy balance: 
dn d£^fc,g , dEth^g dEth^sh fdEk^sh , dEk^RT\ dEg /Air-x 



where 



k,sh 



dt 

1 -^N 



\j:{^^Lu'+2^i,J^^ (AI6) 



dt 



2 



dEth,g _ 



E 



dt ti 7(7 - 1) 



-y2 



PNM^2 I PlM^2 I 7(7 - 1) Phalo ^; 

'-^NM "I '-^/Af "1 7^ '-'halo 

Pa Pa ^ Pa . 



dt 



dEh.RT 1 



Nrt 



(A17) 



(A18) 



-Y^PLUg (A19) 



7^ E Ai^Tw', (A20) 



dt 2 

( |A15| ) gives a good estimate of the radiative shell luminosity if one sets 

dEu^ = 0. (A21) 
dt ^ ^ 

However, in the thin layer approximation there is no satisfactory procedure for the 
calculation of an adiabatic shell thermal energy. Therefore for quasi-adiabatic shell sectors 
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we use: 

Na 

Li = ^j: nl^MTsk,a)ARdJ:^, (A22) 
1=1 

where ^ is the galactic gas metaUicity, At(T) is coohng function, nsh,a and Tsh,a are the mean 
atomic particle number density and temperature within the adiabatic part of the shell, 
AR = 0.14i?^ is the quasi-adiabatic shell thickness, dS"^ is the surface area of the adiabatic 
Lagrangian mesh, and Na is the number of the " adiabatic" Lagrangian meshpoints. 

= JI-MdS^ ^^^^^ 

1^ "'sh,a 

The total shell luminosity then follows from the expression 

L,, = + Lf„ (A25) 
whereas the bolometric bubble luminosity is 

Lboi = Lin + Lsh, (A26) 
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Appendix B: The inner bubble structure 

Our approximate description of the bubble structure is based on the similarity solution 
of Weaver et al. (1977). We assume that the density and temperature distribution can be 
approximated by: 

n = nc{l-x,y\ (Bl) 

T = T,{l-x,)\ (B2) 

where Xi = r/Ri is the dimensionless distance from the cavity center to the particular 
Lagrangian element i, and A = 2/5. To determine ric (and Tc) we assume that the total 
mass interior to superbubble at all times is: 

Mtot = MsN + Me,, (B3) 

where 

MsN= / M^Tvdt, (B4) 



Me. = f Mev^t. (B5) 
Jo 

The value of Msn is defined by (|20|) , and the shell thermal evaporation is (Mac Low & 
McCray, 1988): 

r ^E^cTj'/'^Y.f^^^, radiative phase, 
Me, = j ' (B6) 

I 0, adiabatic phase, 

where k is the Boltzmann's constant, C = 6 x 10"'' erg s~^ cm"^ K"^/^ is the constant 
thermal conductivity coefficient for a fully ionized hydrogen plasma (Cowie & McKee 1977). 

Two approximations are then used. We assume that a portion of the total mass (Mi) 
is enclosed within the radius of the contact discontinuity {Rcd = £iRs) 

p2n rn fRcD . „ 

Ml = i^inu, / d(f) / sinM^ / (1 - r/i?,)-Vdr = 3/i,„ne/^i(A)fi, (B7) 
Jo Jo Jo 
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/„^i(A) 
1 



(1 - a;)"Vdx = 
1 — A 



3- A 



l-(l-£i) 



3-A 



2- A 



-(l-^l)^-' 



(B8) 



Ml — Mtot during the adiabatic phase but as sectors of the shock become radiative, Mtot 
fills a larger volume. In such case, the fraction of Mtot that now reaches up to Rg, assuming 
a collapsed thin shell with radius equal to that of the shock, is simply 



M2 =£2(M-Mi), 



(B9) 



M = 3^,nnJm{X)^, (BIO) 
UA)^/'(l-.)-Vd.^^ + 3l^-jl^. (BID 

The mass distribution of the (adiabatic-radiative) remnant interior can be defined as 

Ml + M2 = e2M + (1 - £2)Mi = Mtot. (B12) 

Equations (B3 - B12) define the central superbubble density and temperature: 

1 Mtot 

'''^"3/.,„[£2/™(A) + (l-£2)/n^l(A)] Q ' ^ 
rjj 3/ii„[£2-^m(A) + (1 — £2)Iml 

^ k mj:;- ^^''^ 

6.0.1. X-ray luminosity 

The total bubble X-ray luminosity includes the X-ray emission from the hot bubble 
interior and the X-ray emission of the swept-up interstellar gas. We conventionally divide 
the X-ray emission from the bubble interior into an emission L^i, from the matter within 
the radius Rcd, and that arising from the gas between the "adiabatic CD" surface and the 
X-ray cutt-off surface, Lx2- The first part is defined by the equation 

^.1 = r d0 r sin ede r""' inl{r)Kr''dr = ZiX-W,SlIxi (B15) 

JO JO Jo 
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where 



„ ~ 1-3A 

/ A.(T) (- 



Tc Jtc 



dT 



The second part we approximate by the value 

Lx2 = {Lxr — Lxi)e2, 



(B16) 



(BIT) 



where L^.^ is the X-ray emission from the interior of the fully radiative bubble (see Silich et 
al. 1996): 

rRs 



1-2-K 




/ 


t) r sin Ode [ 


h 


Jo Jo 



J- C ^^x,cut ^ C 



dT. 



(B18) 
(B19) 



Here Ten is the temperature at the "adiabatic CD" surface, and T^^cut ~ 5 x lO^i^ is the 
X-ray cut-off temperature (see Chu & Mac Low, 1990; also SBHL). The X-ray luminosity 
of the hot bubble interior is then 



+ (1 ~ £2)Lxl- 



(B20) 



The X-ray luminosity from the shell includes both the X-ray emission from the quasi- 
adiabatic sectors of the shell L° and the X-ray emission from the outer layers of the 
radiative shell L^g^- We approximate L'^ ^^ by the expression for the homogeneous shell 
luminosity 

Lth = e E <aA.(T,,,J Ai?dS^, (B21) 

and use the expression 



-R 
-'x,sh 



- x,cuti 
2 



(B22) 



0, T, < Z 

X^Nr ( TjnaA ^ hx(Ts) tR T -> T 

\^ 2^i = l\ Ts ) \tot{Tmax)^b,sh^ J^S^J-X,CUt, 

to estimate the X-ray luminosity from the outer shock front. Here the mean shell density 
and temperature are defined by formulae ( |A23| ) - ( [A24| ) , dS"^ is the surface area of the 
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adiabatic Lagrangian mesh, Tj^ax ~ lO^X is the temperature at the peak of the coohng 
function Afot- The shell X-ray emission then equals 

Lx,out = L^sh + L^^sh (B23) 
and the total bubble X-ray luminosity reads as 

Lx — Lx^in + Lx^out- (B24) 



6.0.2. Phase transitions 

During the evolution in a disk-halo system, superbubbles undergo several transitions 
(from quasi-adiabatic to radiative and viceversa). The time for phase transition can be 
estimated by comparing the characteristic cooling time 

Tcooi = 7, — 7 — X— • B25 

with the dynamic time (t). Here Ua is the atomic number density and the ratio of //a/A*m 
takes into account gas ionization. We assume the mean shell density and temperature (see 
Appendix A) to estimate the time for the adiabatic - radiative transition, and use Gaetz 
and Salpeter (1983) cooling function up to 4 x IQ^K. The T^/^ growth of the cooling rate 
for higher temperatures is assumed. The sudden loss of the thermal energy during the 
adiabatic to radiative transition (see Chevaher 1974; and Falle 1981), is taken into account 
and is assumed to be in a direct proportion to the surface of the rapidly cooling part of the 
shell: 

AEth = jEth. (B26) 

Thus, if the adiabatic-radiative transition occurs for a spherical shell, the remnant looses 
half it's current thermal energy. 
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6.0.3. Cusps formation and grid rearrangement 

As the shell reaches several characteristic scales H^, the acceleration phase begins, and 
the shell experiences Rayleigh- Taylor instabilities. Once the accelerated sector of the shell 
is disrupted, a secondary shell forms in the low density halo. Soon the bottom parts of the 
"cap" begin to overlap and intersect the vertically elongated main body of the remnant and 
a ring-like cusp forms. Some of the Lagrangian particles then move into the bubble interior 
what makes the later calculations impossible. To avoid this difficulty we have introduced 
a simple geometric criterion. The cylindrical radius of any Lagrangian element at any 
Lagrangian "^" -slice (with the exception of top, bottom, and those, which have maximum 
radii) should be either greater than the cylindrical radius of the previous Lagrangian mesh 
element and smaller than the radius of the next one, or it should be smaller and greater 
than the radii of the previous and following Lagrangian mesh particles. If this is not the 
case, we assume cusp formation and the mass of the corresponding particle goes into the 
remnant interior and is not considered any further. 

We also exclude from future consideration those Lagrangian zones which evolve and 
interpenetrate the galaxy midplane. To maintain a sufficient space resolution we regrid the 
shell, if the total number of Lagrangian zones, which are parallel to the galaxy plane falls 
below the initial value N^. 
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Figure captions 

1) The ISM of dwarf galaxies, a-d show the density distribution for models A120, AlOO, 
A80, and A60. For comparison purpose e shows model B60 and f shows the distribution 
obtained for model C50 (see Table 1). In all frames the outer contour was drawn at the 
edge of the galaxy. 

2) The ISM mass distribution. The distribution of the 10^ Mq and 10^ Mq ISM for 
models A and C (see Table 1), respectively. 

3) The time evolution. Evolutionary sequence showing the shape and development of 
the superbubble for the cases AlOO and B60 (see Table 1). The various contours are drawn 
at t = 14 xlO^ yr, 26 xlO^ yr and 50 xlO^ yr. 

4) The time evolution. Panel a displays the amount of ISM swept by the superbubble 
calculated for the case AlOO (see Table 1), as a function of time, while b and c display the 
matter deposited into the superbubble interior by thermal evaporation from the shell of 
swept up matter and that deposited by the sequential SNe. d and e present the calculated 
bolometric and X-ray luminosity of the superbubble. The latter one shows the contribution 
from the swept up shell being much smaller than that produced by the remnant interior 
which accounts for most of the total X-ray luminosity. 

5) The fate of processed matter. Panels a and b compare the maximum expansion 
speed, measured along the symmetry axis (solid lines), with the escape velocity as a 
function of distance to the centre of the galaxy (dashed lines), for all A and B cases. The 
dotted lines mark the moment when the fastest expansion speed drops below the assumed 
halo velocity dispersion. 

6) The fate of processed matter. The ratio of the pressure within the superbubble (for 
models A) to that at the edge of its host galaxy, as a function of the superbubble size. The 
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sudden change in slope is , in all cases, due to the sudden cut in energy injection at tburst — 
40 xlO^ yr. The dotted line indicates the run of the super bubble interior pressure for the 
model A60 normalized to the ISM pressure at the shell stall radius. 

7) The critical luminosity. Panels a and b show the critical luminosity values (derived 
from eq 26; solid lines) required for a superbubble evolving into a Mjsm = 10^ M0 (a) or 
MjsM = 10^ M0 (b) to reach the edge of its host galaxy with a velocity similar to the 
escape velocity measured at the edge of the system. Larger luminosities would produce 
faster remnants and thus their processed matter would be lost from the galaxy. Panel (a) 
also shows the luminosity values obtained from the same eq 26 (dashed line) for all B cases; 
i.e. the cases for which the dark matter component is assumed to be distributed within a 
radius of 5 kpc, instead of 10 kpc. The dotted line (in a and b) indicates the luminosity 
values derived from eq 28, regarded as a better approximation when the evolutionary time 
exceeds the tburst = 40 x 10^ yr. The diamond and triangular symbols indicate the critical 
luminosity values derived from our numerical simulations, after a number of iterations. For 
each of these, the indicated ages mark the evolutionary time required for the superbubble 
to reach the edge of the system. Those exceeding tburst should then be compared with the 
solution represented with the dotted line. 
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Table 1 

Model parameters 





Mtot 


Mgas 


Rb 


Rg 


Chalo 


no 


^burst 


L 


Shell evap. 


N 


lO^Mo 




kpc 


kpc 


km/s~^ 




10^^ ergs 


lO-iOerg s-i 




A120 


10.0 


1.0 


10 


3.9 


120 


20.3 


10 


7.9 


yes 


AlOO 


10.0 


1.0 


10 


8.5 


100 


20.2 


10 


7.9 


yes 


A80 


10.0 


1.0 


10 


13.6 


80 


20.2 


10 


7.9 


yes 


A60 


10.0 


1.0 


10 


24.0 


60 


20.2 


5 


3.95 


yes 


A* 100 


10.0 


1.0 


10 


8.5 


100 


20.2 


10 


7.9 


no 


B160 


10.0 


1.0 


5 


2.8 


160 


20.6 


10 


7.9 


yes 


B120 


10.0 


1.0 


5 


6.0 


120 


20.3 


10 


7.9 


yes 


BlOO 


10.0 


1.0 


5 


8.7 


100 


20.3 


10 


7.9 


yes 


B60 


10.0 


1.0 


5 


24.0 


60 


21.6 


5 


3.95 


yes 


C50 


1.0 


0.1 


5 


3.0 


50 


20.05 


1.07 


0.85 


yes 


C40 


1.0 


0.1 


5 


5.4 


40 


20.20 


0.57 


0.45 


yes 


( 


1.0 


0.1 


.) 


9.7 


:-!() 


20.20 


0. 12 


0.:i:-! 




D80 


5.0 


1.0 


10 


5.7 


80 


20.26 


10 


7.9 


yes 


E80 


5.0 


1.0 


5 


6.8 


80 


20.36 


10 


7.9 


yes 



Model AlOO _ 




R(kpc) 



15 




R(kpc) 



30 



Model A60 




10 20 50 
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